Metagenomic analysis of individual mosquito viromes reveals the geographical patterns and drivers of viral diversity.

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Here, using a meta-transcriptomic approach, we determined the viromes of 2,438 individual mosquitoes (81 species), spanning ~4,000 km along latitudes and longitudes in China. From these data we identified 393 viral species associated with mosquitoes, including 7 (putative) species of arthropod-borne viruses (that is, arboviruses). We identified potential mosquito species and geographic hotspots of viral diversity and arbovirus occurrence, and demonstrated that the composition of individual mosquito viromes was strongly associated with host phylogeny. Our data revealed a large number of viruses shared among mosquito species or genera, enhancing our understanding of the host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, perhaps reflecting long-distance mosquito dispersal. Together, these results greatly expand the known mosquito virome, linked viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the biogeography and diversity of viruses in insect vectors.

Total Infectome Characterization of Respiratory Infections during the 2022-23 COVID-19 Outbreak in China Revealed Extensive Coinfections with Links to SARS-CoV-2 Status, Age, and Disease Severity.

Between 7 December 2022 and 28 February 2023, China experienced a new wave of COVID-19 that swept across the entire country and resulted in an increasing amount of respiratory infections and hospitalizations. The purpose of this study is to reveal the intensity and composition of coinfecting microbial agents. In total, 196 inpatients were recruited from The Third People's Hospital of Shenzhen, and 169 respiratory and 73 blood samples were collected for metagenomic next-generation sequencing. The total "Infectome" was characterized and compared across different groups defined by the SARS-CoV-2 detection status, age groups, and severity of disease. Our results revealed a total of 22 species of pathogenic microbes (4 viruses, 13 bacteria, and 5 fungi), and more were discovered in the respiratory tract than in blood. The diversity of the total infectome was highly distinguished between respiratory and blood samples, and it was generally higher in patients that were SARS-CoV-2-positive, older in age, and with more severe disease. At the individual pathogen level, HSV-1 seemed to be the major contributor to these differences observed in the overall comparisons. Collectively, this study reveals the highly complex respiratory infectome and high-intensity coinfection in patients admitted to the hospital during the period of the 2023 COVID-19 pandemic in China.

The clinical outcome of COVID-19 is strongly associated with microbiome dynamics in the upper respiratory tract.

OBJECTIVES: The respiratory tract is the portal of entry for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although a variety of respiratory pathogens other than SARS-CoV-2 have been associated with severe cases of COVID-19 disease, the dynamics of the upper respiratory microbiota during disease the course of disease, and how they impact disease manifestation, remain uncertain. METHODS: We collected 349 longitudinal upper respiratory samples from a cohort of 65 COVID-19 patients (cohort 1), 28 samples from 28 recovered COVID-19 patients (cohort 2), and 59 samples from 59 healthy controls (cohort 3). All COVID-19 patients originated from the earliest stage of the epidemic in Wuhan. Based on a modified clinical scale, the disease course was divided into five clinical disease phases (pseudotimes): "Healthy" (pseudotime 0), "Incremental" (pseudotime 1), "Critical" (pseudotime 2), "Complicated" (pseudotime 3), "Convalescent" (pseudotime 4), and "Long-term follow-up" (pseudotime 5). Using meta-transcriptomics, we investigated the features and dynamics of transcriptionally active microbes in the upper respiratory tract (URT) over the course of COVID-19 disease, as well as its association with disease progression and clinical outcomes. RESULTS: Our results revealed that the URT microbiome exhibits substantial heterogeneity during disease course. Two clusters of microbial communities characterized by low alpha diversity and enrichment for multiple pathogens or potential pathobionts (including Acinetobacter and Candida) were associated with disease progression and a worse clinical outcome. We also identified a series of microbial indicators that classified disease progression into more severe stages. Longitudinal analysis revealed that although the microbiome exhibited complex and changing patterns during COVID-19, a restoration of URT microbiomes from early dysbiosis toward more diverse status in later disease stages was observed in most patients. In addition, a group of potential pathobionts were strongly associated with the concentration of inflammatory indicators and mortality. CONCLUSION: This study revealed strong links between URT microbiome dynamics and disease progression and clinical outcomes in COVID-19, implying that the treatment of severe disease should consider the full spectrum of microbial pathogens present.

ICTV Virus Taxonomy Profile: Jingchuvirales 2023.

Jingchuvirales is an order of negative-sense RNA viruses with genomes of 9.1-15.3 kb that have been associated with arachnids, barnacles, crustaceans, insects, fish and reptiles in Africa, Asia, Australia, Europe, North America and South America. The jingchuviral genome has two to four open reading frames (ORFs) that encode a glycoprotein (GP), a nucleoprotein (NP), a large (L) protein containing an RNA-directed RNA polymerase (RdRP) domain, and/or proteins of unknown function. Viruses in the order are only known from their genome sequences. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the order Jingchuvirales and on the families Aliusviridae, Chuviridae, Crepuscuviridae, Myriaviridae and Natareviridae, which are available at ictv.global/report/jingchuvirales, ictv.global/report/aliusviridae, ictv.global/report/chuviridae, ictv.global/report/crepuscuviridae, ictv.global/report/myriaviridae and ictv.global/report/natareviridae, respectively.

Refining PRRSV-2 genetic classification based on global ORF5 sequences and investigation of their geographic distributions and temporal changes.

IMPORTANCE: In this study, comprehensive analysis of 82,237 global porcine reproductive and respiratory syndrome virus type 2 (PRRSV-2) open reading frame 5 sequences spanning from 1989 to 2021 refined PRRSV-2 genetic classification system, which defines 11 lineages and 21 sublineages and provides flexibility for growth if additional lineages, sublineages, or more granular classifications are needed in the future. Geographic distribution and temporal changes of PRRSV-2 were investigated in detail. This is a thorough study describing the molecular epidemiology of global PRRSV-2. In addition, the reference sequences based on the refined genetic classification system are made available to the public for future epidemiological and diagnostic applications worldwide. The data from this study will facilitate global standardization and application of PRRSV-2 genetic classification.

Origin, genomic diversity and evolution of African swine fever virus in East Asia.

Since 2018, the outbreaks of genotype II African swine fever virus (ASFV) in China and several eastern Asian countries have caused a huge impact on the local swine industry, resulting in huge economic losses. However, little is known about the origin, genomic diversity, evolutionary features, and epidemiological history of the genotype II ASFV. Here, 14 high-quality complete genomes of ASFVs were generated via sequencing of samples collected from China over the course of 3 years, followed by phylogenetic and phylodynamic analyses. The strains identified were relatively homogeneous, with a total of 52 SNPs and 11 indels compared with the prototype strain HLJ/2018, among which there were four exceptionally large deletions (620-18,023 nt). Evolutionary analyses revealed that ASFV strains distributed in eastern Asia formed a monophyly and a 'star-like' structure centered around the prototype strain, suggesting a single origin. Additionally, phylogenetic network analysis and ancestral reconstruction of geographic state indicated that genotype II ASFV strains in eastern Asia likely originated from Western Europe. Overall, these results contribute to the understanding of the history and current status of genotype II ASFV strains in eastern Asian, which could be of considerable importance in disease control and prevention.

Viral community composition of hypersaline lakes.

Despite their widespread distribution and remarkable antiquity no RNA viruses definitively associated with the domain Archaea have been identified. In contrast, 17 families of DNA viruses are known to infect archaea. In an attempt to uncover more of the elusive archaeal virosphere, we investigated the metatranscriptomes of hypersaline lakes that are a rich source of archaea. We sequenced RNA extracted from water filter samples of Lake Tyrrell (Victoria, Australia) and cultures seeded from four lakes in Antarctica. To identify highly divergent viruses in these data, we employed a variety of search tools, including Hidden Markov models (HMMs) and position-specific scoring matrices (PSSMs). From this, we identified 12 highly divergent, RNA virus-like candidate sequences from the virus phyla Artverviricota, Duplornaviricota, Kitrinoviricota, Negarnaviricota, and Pisuviricota, including those with similarity to the RNA-dependent RNA polymerase (RdRp). An additional analysis with an artificial intelligence (AI)-based approach that utilises both sequence and structural information identified seven putative and highly divergent RdRp sequences of uncertain phylogenetic position. A sequence matching the Pisuviricota from Deep Lake in Antarctica had the strongest RNA virus signal. Analyses of the dinucleotide representation of the virus-like candidates in comparison to that of potential host species were in some cases compatible with an association to archaeal or bacterial hosts. Notably, however, the use of archaeal CRISPR spacers as a BLAST database failed to detect any RNA viruses. We also described DNA viruses from the families Pleolipoviridae, Sphaerolipoviridae, Halspiviridae, and the class Caudoviricetes. Although we were unable to provide definitive evidence the existence of an RNA virus of archaea in these hypersaline lakes, this study lays the foundations for further investigations of highly divergent RNA viruses in natural environments.

Metagenomic analysis of individual mosquitos reveals the ecology of insect viruses.

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Using a meta-transcriptomic approach, we analysed the virome of 2,438 individual mosquitos (79 species), spanning ~4000 km along latitudes and longitudes in China. From these data we identified 393 core viral species associated with mosquitos, including seven (putative) arbovirus species. We identified potential species and geographic hotspots of viral richness and arbovirus occurrence, and demonstrated that host phylogeny had a strong impact on the composition of individual mosquito viromes. Our data revealed a large number of viruses shared among mosquito species or genera, expanding our knowledge of host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, possibly facilitated by long-distance mosquito migrations. Together, our results greatly expand the known mosquito virome, linked the viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the ecology of viruses of insect vectors.

Association of vaccination, international travel, public health and social measures with lineage dynamics of SARS-CoV-2.

Continually emerging SARS-CoV-2 variants of concern that can evade immune defenses are driving recurrent epidemic waves of COVID-19 globally. However, the impact of measures to contain the virus and their effect on lineage diversity dynamics are poorly understood. Here, we jointly analyzed international travel, public health and social measures (PHSM), COVID-19 vaccine rollout, SARS-CoV-2 lineage diversity, and the case growth rate (GR) from March 2020 to September 2022 across 63 countries. We showed that despite worldwide vaccine rollout, PHSM are effective in mitigating epidemic waves and lineage diversity. An increase of 10,000 monthly travelers in a single country-to-country route between endemic countries corresponds to a 5.5% (95% CI: 2.9 to 8.2%) rise in local lineage diversity. After accounting for PHSM, natural immunity from previous infections, and waning immunity, we discovered a negative association between the GR of cases and adjusted vaccine coverage (AVC). We also observed a complex relationship between lineage diversity and vaccine rollout. Specifically, we found a significant negative association between lineage diversity and AVC at both low and high levels but not significant at the medium level. Our study deepens the understanding of population immunity and lineage dynamics for future pandemic preparedness and responsiveness.

A metatranscriptomic analysis of geothermal hot springs reveals diverse RNA viruses including the phylum Lenarviricota.

Little is known about the diversity of RNA viruses in geothermal systems. We generated total RNA sequencing data from two hot springs in Kuirau Park, Rotorua, New Zealand. In one data set, from a 71.8 °C pool, we observed a microbial community that was 98.5% archaea. The second data set, representing a cooler 36.8 °C geothermal hot spring, had a more diverse microbial profile: 58% bacteria, 34.5% eukaryotes and 7.5% archaea. Within this latter pool, we detected sequences likely representing 23 RNA viruses from the families Astroviridae, Tombusviridae, Polycipiviridae, Discistroviridae, Partitiviridae, and Mitoviridae, as well as from unclassified clades of the orders Tolivirales, Picornavirales, and Ghabrivirales. Most viruses had uncertain host associations. Of particular note, we identified four novel RNA viruses from the phylum Lenarviricota, commonly associated with bacteria and fungi, that occupied a divergent phylogenetic position within unclassified clades and may represent an ancient order-level taxon of unknown host association.

Individual bat virome analysis reveals co-infection and spillover among bats and virus zoonotic potential.

Bats are reservoir hosts for many zoonotic viruses. Despite this, relatively little is known about the diversity and abundance of viruses within individual bats, and hence the frequency of virus co-infection and spillover among them. We characterize the mammal-associated viruses in 149 individual bats sampled from Yunnan province, China, using an unbiased meta-transcriptomics approach. This reveals a high frequency of virus co-infection (simultaneous infection of bat individuals by multiple viral species) and spillover among the animals studied, which may in turn facilitate virus recombination and reassortment. Of note, we identify five viral species that are likely to be pathogenic to humans or livestock, based on phylogenetic relatedness to known pathogens or in vitro receptor binding assays. This includes a novel recombinant SARS-like coronavirus that is closely related to both SARS-CoV and SARS-CoV-2. In vitro assays indicate that this recombinant virus can utilize the human ACE2 receptor such that it is likely to be of increased emergence risk. Our study highlights the common occurrence of co-infection and spillover of bat viruses and their implications for virus emergence.

Hierarchical metacommunity structure of fungal endophytes.

The ecological and evolutionary processes shaping community structure and functions of microbial symbionts are known to be scale-dependent. Nonetheless, understanding how the relative importance of these processes changes across spatial scales, and deciphering the hierarchical metacommunity structure of fungal endophytes has proven challenging. We investigated metacommunities of endophytic fungi within leaves of an invasive plant (Alternanthera philoxeroides) across wide latitudinal transects both in its native (Argentina) and introduced (China) ranges to test whether metacommunities of fungal endophytes were structured by different drivers at different spatial scales. We found Clementsian structures with seven discrete compartments (distinctive groups of fungal species with coincident distribution ranges), which coincided with the distribution of major watersheds. Metacommunity compartments were explicitly demarcated at three spatial scales, that is, the between-continent, between-compartment, and within-compartment scales. At larger spatial scales, local environmental conditions (climate, soil, and host plant traits) were replaced by other geographical factors as principal determinants of metacommunity structure of fungal endophytes and community diversity-function relationships. Our results reveal novel insights into the scale dependency of diversity and functions of fungal endophytes, which are likely similar for plant symbionts. These findings can potentially improve our understanding of the global patterns of fungal diversity.

Nelson Bay Reovirus Isolated from Bats and Blood-Sucking Arthropods Collected in Yunnan Province, China.

Nelson Bay reovirus (NBV) is an emerging zoonotic virus that can cause acute respiratory disease in humans. These viruses are mainly discovered in Oceania, Africa, and Asia, and bats have been identified as their main animal reservoir. However, despite recent expansion of diversity for NBVs, the transmission dynamics and evolutionary history of NBVs are still unclear. This study successfully isolated two NBV strains (MLBC1302 and MLBC1313) from blood-sucking bat fly specimens (Eucampsipoda sundaica) and one (WDBP1716) from the spleen specimen of a fruit bat (Rousettus leschenaultii), which were collected at the China-Myanmar border area of Yunnan Province. Syncytia cytopathic effects (CPE) were observed in BHK-21 and Vero E6 cells infected with the three strains at 48 h postinfection. Electron micrographs of ultrathin sections showed numerous spherical virions with a diameter of approximately 70 nm in the cytoplasm of infected cells. The complete genome nucleotide sequence of the viruses was determined by metatranscriptomic sequencing of infected cells. Phylogenetic analysis demonstrated that the novel strains were closely related to Cangyuan orthoreovirus, Melaka orthoreovirus, and human-infecting Pteropine orthoreovirus HK23629/07. Simplot analysis revealed the strains originated from complex genomic reassortment among different NBVs, suggesting the viruses experienced a high reassortment rate. In addition, strains successfully isolated from bat flies also implied that blood-sucking arthropods might serve as potential transmission vectors. IMPORTANCE Bats are the reservoir of many viral pathogens with strong pathogenicity, including NBVs. Nevertheless, it is unclear whether arthropod vectors are involved in transmitting NBVs. In this study, we successfully isolated two NBV strains from bat flies collected from the body surface of bats, which implies that they may be vectors for virus transmission between bats. While the potential threat to humans remains to be determined, evolutionary analyses involving different segments revealed that the novel strains had complex reassortment histories, with S1, S2, and M1 segments highly similar to human pathogens. Further experiments are required to determine whether more NBVs are vectored by bat flies, their potential threat to humans, and transmission dynamics.

Genetic Diversity and Evolution of Viruses Infecting Felis catus: A Global Perspective.

Cats harbor many important viral pathogens, and the knowledge of their diversity has been greatly expanded thanks to increasingly popular molecular sequencing techniques. While the diversity is mostly described in numerous regionally defined studies, there lacks a global overview of the diversity for the majority of cat viruses, and therefore our understanding of the evolution and epidemiology of these viruses was generally inadequate. In this study, we analyzed 12,377 genetic sequences from 25 cat virus species and conducted comprehensive phylodynamic analyses. It revealed, for the first time, the global diversity for all cat viruses known to date, taking into account highly virulent strains and vaccine strains. From there, we further characterized and compared the geographic expansion patterns, temporal dynamics and recombination frequencies of these viruses. While respiratory pathogens such as feline calicivirus showed some degree of geographical panmixes, the other viral species are more geographically defined. Furthermore, recombination rates were much higher in feline parvovirus, feline coronavirus, feline calicivirus and feline foamy virus than the other feline virus species. Collectively, our findings deepen the understanding of the evolutionary and epidemiological features of cat viruses, which in turn provide important insight into the prevention and control of cat pathogens.

Spike substitution T813S increases Sarbecovirus fusogenicity by enhancing the usage of TMPRSS2.

SARS-CoV Spike (S) protein shares considerable homology with SARS-CoV-2 S, especially in the conserved S2 subunit (S2). S protein mediates coronavirus receptor binding and membrane fusion, and the latter activity can greatly influence coronavirus infection. We observed that SARS-CoV S is less effective in inducing membrane fusion compared with SARS-CoV-2 S. We identify that S813T mutation is sufficient in S2 interfering with the cleavage of SARS-CoV-2 S by TMPRSS2, reducing spike fusogenicity and pseudoparticle entry. Conversely, the mutation of T813S in SARS-CoV S increased fusion ability and viral replication. Our data suggested that residue 813 in the S was critical for the proteolytic activation, and the change from threonine to serine at 813 position might be an evolutionary feature adopted by SARS-2-related viruses. This finding deepened the understanding of Spike fusogenicity and could provide a new perspective for exploring Sarbecovirus' evolution.

Dysbiosis of oropharyngeal microbiome and antibiotic resistance in hospitalized COVID-19 patients.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic is ongoing and multiple studies have elucidated its pathogenesis, however, the related- microbiome imbalance caused by SARS-CoV-2 is still not clear. In this study, we have comprehensively compared the microbiome composition and associated function alterations in the oropharyngeal swabs of healthy controls and coronavirus disease 2019 (COVID-19) patients with moderate or severe symptoms by metatranscriptomic sequencing. We did observe a reduced microbiome alpha-diversity but significant enrichment of opportunistic microorganisms in patients with COVID-19 compared with healthy controls, and the microbial homeostasis was rebuilt following the recovery of COVID-19 patients. Correspondingly, less functional genes in multiple biological processes and weakened metabolic pathways such as carbohydrate metabolism, energy metabolism were also observed in COVID-19 patients. We only found higher relative abundance of limited genera such as Lachnoanaerobaculum between severe patients and moderate patients while no worthy-noting microbiome diversity and function alteration were observed. Finally, we noticed that the co-occurrence of antibiotic resistance and virulence was closely related to the microbiome alteration caused by SRAS-CoV-2. Overall, our findings demonstrate that microbial dysbiosis may enhance the pathogenesis of SARS-CoV-2 and the antibiotics treatment should be critically considered.

Zoonotic risk assessment among farmed mammals.

Farmed mammals may act as hosts for zoonotic viruses that can cause disease outbreaks in humans. This SnapShot shows which farmed mammals, and to what extent, are of particular risk of harboring and spreading viruses from viral families that are commonly associated with zoonotic disease. It also discusses genome surveillance methods and biosafety measures. To view this SnapShot, open or download the PDF.

Emergence of an ancient and pathogenic mammarenavirus.

Emerging zoonoses of wildlife origin caused by previously unknown agents are one of the most important challenges for human health. The Qinghai-Tibet Plateau represents a unique ecological niche with diverse wildlife that harbours several human pathogens and numerous previously uncharacterized pathogens. In this study, we identified and characterized a novel arenavirus (namely, plateau pika virus, PPV) from plateau pikas (Ochotona curzoniae) on the Qinghai-Tibet Plateau by virome analysis. Isolated PPV strains could replicate in several mammalian cells. We further investigated PPV pathogenesis using animal models. PPV administered via an intraventricular route caused trembling and sudden death in IFNαßR-/- mice, and pathological inflammatory lesions in brain tissue were observed. According to a retrospective serological survey in the geographical region where PPV was isolated, PPV-specific IgG antibodies were detected in 8 (2.4%) of 335 outpatients with available sera. Phylogenetic analyses revealed that this virus was clearly separated from previously reported New and Old World mammarenaviruses. Under the co-speciation framework, the estimated divergence time of PPV was 77-88 million years ago (MYA), earlier than that of OW and NW mammarenaviruses (26-34 MYA).

VirBot: an RNA viral contig detector for metagenomic data.

SUMMARY: Without relying on cultivation, metagenomic sequencing greatly accelerated the novel RNA virus detection. However, it is not trivial to accurately identify RNA viral contigs from a mixture of species. The low content of RNA viruses in metagenomic data requires a highly specific detector, while new RNA viruses can exhibit high genetic diversity, posing a challenge for alignment-based tools. In this work, we developed VirBot, a simple yet effective RNA virus identification tool based on the protein families and the corresponding adaptive score cutoffs. We benchmarked it with seven popular tools for virus identification on both simulated and real sequencing data. VirBot shows its high specificity in metagenomic datasets and superior sensitivity in detecting novel RNA viruses. AVAILABILITY AND IMPLEMENTATION: https://github.com/GreyGuoweiChen/RNA_virus_detector. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.